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We study the hadron-quark phase transition with the finite size effects at finite temperature. For 
the hadron phase, we adopt a realistic equation of state in the framework of the Brueckner-Hartree- 
Fock theory including hyperons. The properties of the mixed phase are clarified by considering the 
finite size effects under the Gibbs conditions. We find that the equation of state becomes softer 
than that at zero-temperature for some density region. We also find that the equation of state gets 
closer to that given by the Maxwell construction. Moreover, the number of hyperons is suppressed 
by the presence of quarks. These are characteristic features of the hadron-quark mixed phase, and 
should be important for many astrophysical phenomena such as mergers of binary neutron stars. 

PACS numbers: 97.60.Jd, 12.39.Ba 



I. INTRODUCTION 

Nowadays the effects of quark matter on various astro- 
physical phenomena have been studied extensively. For 
example, the cooling of compact stars have been studied 
in Ref. [H, [2, 0, H|. Other examples include the effects 
of quark matter on gravitational wave radiation [H, H J7 1 , 
neutrino emissions [g, H, flo| , rotational frequencies 11 1 , 
and the energy release during the collapse from neutron 
stars to quark stars [H|,[lH, etc.. In particular, the mech- 
anisms of supernovae and gamma-ray bursts have not 
been clearly understood; the QCD phase transition may 
take place in such phenomena and take critical role [14| . 

However, there are left many uncertainties for the 
hadron-quark phase transition, e.g. the equation of state 
(EOS) of quark matter or dcconffncmcnt mechanism. As- 
suming the quark deconfinement transition to be of first 
order, we find it causes a thermodynamical instability 
and the mixed phase appears around the critical den- 
sity. Since there are two conserved quantities, baryon 
number and electric charge, the phase equilibrium in the 
mixed phase must be carefully treated by applying the 
Gibbs conditions [Tj|, instead of the Maxwell construc- 
tion. A simple treatment of the mixed phase may be 
the bulk Gibbs calculation, where phase equilibrium of 
two bulk matter is considered without electromagnetic 
interaction and surface tension. Generally the properties 
of the mixed phase should strongly depend on electro- 
magnetic interaction and surface tension, and these ef- 
fects, sometimes called "the finite-size effects" , lead to 
the non-uniform "Pasta" structures. The EOS of the 
mixed phase becomes similar to the one under the bulk 
Gibbs calculation for weak surface tension, and to the 
one given by the Maxwell construction for strong surface 
tension [l|| [l?], [HI • The charge screening is also impor- 
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tant for their mechanical instability. 

In the previous papers these finite-size effects have 
been properly taken into account to elucidate the prop- 
erties of the pasta structure and demonstrate the impor- 
tance of the charge screening at zero temperature |18j |. 
However, finite temperature comes in many cases such 
as relativistic heavy-ion collisions and astrophysical phe- 
nomena. In this paper, we study the hadron-quark mixed 
phase with the finite-size effects at finite temperatur e by 
extending the previous works of Maruyama et al. [l8| . 
We adopt the Brueckner-Hartree-Fock EOS by Baldo et 
al. for the hadron phase [lji[ . The EOS includes hyperons 
as well as nucleons, interacting with hadronic two-body 
forces and nucleonic three-body forces. For the quark 
phase, we adopt the thermodynamic bag model for sim- 
plicity. We impose the Gibbs conditions on the phase 
equilibrium, taking into account the finite-size effects. 

This paper is organized as follows. In Sec. II, we out- 
line our framework. In Sec. Ill, we present numerical 
results. Sec. IV is devoted to the conclusion and discus- 
sion where we give some astrophysical implications of our 
results. 



II. EQUATION OF STATE 

A. Equation of state for hadron phase at finite 
temperature — Brueckner-Hartree-Fock theory 

Our theoretical framework for the hadron phase 
of matter is the nonrclativistic Brucckncr-Hartrce- 
Fock approach (20| based on the microscopic nucleon- 
nucleon (NN), nucleon-hyperon (NY), and hyperon- 
hyperon (YY) potentials. The Brueckner-Hartree-Fock 
calculation is a reliable and well-controlled theoretical 
approach for the study of dense baryonic matter. De- 
tailed procedure can be found in Rcfs. [H, HH, HH . 

We adopt the Argonne Vis potential Q for NN po- 
tentials, and semi-phenomcnological Urbana UIX nuclc- 
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onic three body forces 1241 and the Nijmegen soft-core 
NSC89 NY potentials j25[. Unfortunately, there are not 
reliable YY potentials now, because no YY scattering 
data exist. Therefore, we neglect YY potentials in this 
paper. Recently, the maximum mass of neutron stars are 
calculated ^ using the NSC97 NY and YY potentials 
[2^ . However, their maximum masses are quite similar. 

In this paper, we discuss the hadron-quark mixed 
phase with finite-size effects at finite temperature [H3 [. 
At finite temperature, there is no microscopic calculation. 
To take into account the effects of finite tem per ature, we 
adopt Frozen Correlations Approximation [28l l29l [30| . 
In this approximation, the correlations at finite temper- 
ature are assumed to be the same with the ones at zero 
temperature. It is found to be good accuracy at finite 
temperature by past studies [H, [29|, [3(| • Accordingly, we 
focus on hadron-quark mixed phase at the finite temper- 
ature and adopt this approximation for the neutrino-free 
case. 

In the following, we briefly describe our framework. 
First we get the chemical potential /ij from the number 
density rii, 
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where Si and fi (p) are the single-particle energy and the 
Fermi-Dirac distribution function, respectively, whereas 
subscript i shows the particle species, i = n,p, A, £~. We 
set each degeneracy factor g = 2, and adopt each mass as 
m n = m. p = 939 MeV, ttia = 1115.7 MeV, and m s - = 
1197.4 MeV. We note that £j includes the interaction 
energy Ui as well as the kinetic energy [28l f3l| . 



e t = yJml+pP + Ui. (3) 
Finally, we get the free-energy density T as following 
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B. Equation of state for quark phase 
— The thermodynamic bag model 

Unfortunately, the current theoretical description of 
quark matter includes many uncertainties, seriously lim- 
iting the predictability of EOS at high baryon density. 
For the time being, we only resort to a phenomenological 
model for the quark matter EOS and try to put the con- 
straint of parameters by a few experimental information. 

In this paper, we adopt the thermodynamic bag model 
to construct EOS for simplicity. The bag model may be 
an oversimplified model for quark matter. It should be 
interesting to compare our results with those given by 
other models; e.g. Nambu-Jona-Lasinio model [HI, [33[ 
or Polyakov-Nambu- Jona-Lasinio model [H], [35[ . 

We get the number density uq, the pressure Pq, and 
the energy density 6q, 
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for uniform quark matter, where f q (p) is the Fermi-Dirac 
distribution function of the quark q (= u,d,s), m q is 
its current mass, and B is the energy density difference 
between the perturbative vacuum and the true vacuum, 
the bag constant. In this paper, we employ B = 100 MeV 
fm -3 , and the degeneracy factor g = 6. We use massless 
u and d quarks and m s = 150 MeV, and ignore any 
anti-particles because quark chemical potential are much 
larger than masses and temperature. 
The baryon number is conserved as 
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where sh is the entropy density calculated from 
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by 



The total pressure for uniform hadron phase is given 



Ph = E' ij71 ' ~ 



(6) 



C. Hadron-quark mixed phase under the Gibbs 
conditions 

To take into account the finite-size effects, we impose 
the Gibbs conditions on the mixed phase [361 ] . which re- 
quire the pressure balance and the equality of the chem- 
ical potentials between two phases besides the thermal 
equilibrium. We employ the Wigner-Seitz approxima- 
tion in which the whole space is divided into equivalent 
cells with given geometrical symmetry, specified by the 
dimensionality d — 3 (droplet or bubble), d = 2 (rod 
or tube), or d = 1 (slab). The structures of tube and 
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bubble are opposite distributions of rod and droplet as 
illustrated in Fig. [1] 
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FIG. 1: Pasta structures in the Wigner-Seitz approximation. 
We assume that one of them appears in the hadron-quark 
mixed phase under the periodic condition. 

The quark and hadron phases are separated in each cell 
with volume Vw ■ a lump made of the quark phase with 
volume Vq is embedded in the hadronic phase with vol- 
ume Vh or vice versa. A sharp boundary is assumed be- 
tween the two phases and the surface energy is taken into 
account in terms of a surface-tension parameter a. The 
surface tension of the hadron-quark interface is poorly 
known, but some theoretical estimates based on the MIT 



bag model for strangelets 37 1 and lattice gauge simula- 
tions at finite temperature 38( suggest a range of a w 10- 
100 MeV/fm 2 . We show results using a = 40 MeV/fm 2 
in this article. 

We use the Thomas-Fermi approximation for the den- 
sity profiles of hadrons and quarks. The Hclmholtz free 
energy for each cell is then given as 



t (r)]+F e + E c + E s 



F= / d 3 rF H [rii(r)}+ f d 3 r T Q [ 
Jv H Jv q 

(13) 

with i = n,p,A,T, , q = u,d,s, Th (J~q) is the 
Helmholtz free energy density for hadron (quark) mat- 
ter, and Es = aS the surface energy with S being the 
hadron-quark interface area. F e is the free energy of the 
electron gas. For simplicity, muons are not included in 
this paper. The value of Eq is the Coulomb interaction 
energy calculated by, 
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where the charge density is given by 
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in Vq with Qi (or Q q ) being the particle charge (Q = 
— e < for the electron). Accordingly, the Coulomb po- 
tential (j>{r) is defined as 
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where cf>o is an arbitrary constant representing the gauge 
degree of freedom. We fix it b y st ipulating the condition, 
(j>{Rw) = 0, as before (39l. |40|. I4l1| . The Poisson equation 
then reads 
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Under the Gibbs conditions, we must consider chemi- 
cal equilibrium at the hadron-quark interface as well as 
inside each phase; 



fl u + fi e = fl d = fis , 

fip + Me = Mn = fJ-A = Mu + ZfJ-d 

MS- + = • 

For a given baryon number density 
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Eqs. p^ri^l) . together with the global charge neutrality 
condition, J Vw d 3 rn^r) = 0, obviously fulfill the require- 
ment by the Gibbs conditions. 



III. NUMERICAL RESULTS 

A. Thermal effects on EOS 

Using above relations, we study the hadron-quark 
mixed phase at finite temperature. Upper panels of 
Fig. [2] indicate the resulting pressures of the hadron- 
quark mixed phase with that of the pure hadron and 
quark phases over the relevant range of baryon den- 
sity at zero temperature (the upper left panel) and fi- 
nite temperature (the upper right panel). Here, we take 
T = 30 McV. The thin curves indicate the pure hadron 
and quark phases, while the thick dots indicate the mixed 
phase in its various geometric realizations by the full cal- 
culation. The transitions between the different geomet- 
rical structures are, by construction, discontinuous and 
a more sophisticated approach would be required for a 
more realistic description of this feature. We also show, 
for comparison, the hadron-quark phase transition by the 
Maxwell construction by the thin solid line. 
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FIG. 2: (Color online) EOS of the mixed phase (thick dots) in comparison with pure hadron and quark phases (thin curves). 
The upper panels show the pressure at zero temperature (left), and T — 30 MeV (right). We also show, for comparison, the 
mixed phase by the Maxwell construction by thin solid line. Lower panels: The free energy per baryon F/A of the mixed phase 
(thick curves) in comparison with pure hadron and quark phases (thin curves). The left panels show the zero temperature case, 
and the right panels the finite temperature case (T — 30 MeV). Different colored segments of the mixed phase are chosen by 
minimizing the energy. 



Compared with the zero temperature case, the mixed 
phase is restricted and EOS gets close to that by the 
Maxwell construction, though we properly apply the 
Gibbs conditions. The restriction of the mixed phase has 
been already demonstrated at T = due to the charge 
screening effect [l8j]. We can see that the geometrical 
structure also becomes unstable due to the thermal ef- 
fects. The left panel of Fig. [3] shows the free energies per 
baryon of the droplet structure at several values of tem- 
perature. The quark volume fraction (R/Rw) 3 is fixed 
to exclude the trivial R— dependence. Here we use, for 
example, the optimal value of (R/Rw) 3 at T = MeV 
in every curve. We normalize them by subtracting the 
free energy at infinite radius, AF = F(R) — F(R — > oo), 
to show the R dependence clearly. The structure of the 
mixed phase is mechanically stable below T ~ 60 MeV, 
but the optimal value of the radius R is shifted to the 
larger value as T increases. This behavior is a signal 
of the mechanical instability and comes from the charge 
screening effect and the thermal effect. 



To elucidate this point more clearly, we show the each 
contribution to AF/A in the right panel of Fig. [3J i.e. 
the Coulomb energy per baron Eq /A, the surface energy 
per baryon Eg /A, and the correlation energy per baryon 
Ecoxt/A. The baryon density and the temperature are 
set as riB = 2no and T = 50 MeV. 

When we treat the Coulomb potential and charge den- 
sities in a self-consistent manner, we can see the charge 
screening effect; it gives rise to the Debye screening mass 
for the Coulomb interaction and induces the rearrange- 
ment of charge densities (see Fig. 4). Then Eq/A expo- 
nentially decreases as R increases for large R. At finite 
temperature the value of Ec decreases due to the hy- 
peron mixing in comparison with the zero temperature 
case, but it is a little effect. On the other hand ther- 
mal effect mainly emerges through the rearrangement of 
charge densities; the kinetic and strong-interaction ener- 
gies give rise to a new R dependence, which is called the 
correlation energy in Ref. [40j , proportional to the differ- 
ence of charge configuration from the uniform one in the 
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absence of the Coulomb interaction. It has been shown 
that E COTT /A exhibits — R^ 1 behavior for large R [io| . 
Its magnitude decreases as T increases due to the re- 
duction of the Coulomb interaction. Thus the optimal 
value of R is shifted to a larger value at finite tempera- 
ture. Note that the extreme case such that R(Rw) — > oo 
corresponds to the Maxwell construction for bulk mat- 
ter. In our formulation, the pasta structures disappear 
at T > 60 McV. 

The pressure of the mixed phase at finite temperature 
is lower than that at zero temperature as shown in the 
upper right panel of Fig. [21 while the pressure of the pure 
hadron or quark phase at finite temperature becomes 
higher than that at zero temperature. When we write the 
pressure P as the sum of the cold part (zero temperature) 
Pcoid and the thermal part Pthermai, P — Pcold ~t~ P thermal 3 
usually we find Pthermai > 0. On the other hand, we find 
-Pthermai < for the mixed phase. Similar results have 
been shown by previous studies, though they did not in- 
clude hyperons nor the finite-size effects [33. H3|- This 
behavior is characteristic for the mixed phase. How does 
such characteristic behavior of pressure appear? 

Since the degrees of freedom for quark matter is larger 
than those for hadron matter, the change of heat in quark 
matter becomes larger than that in hadron matter as 
temperature is increased. Hence the free energy change of 
quark matter is larger than that of hadron matter (com- 
pare the lower panels of Fig. [J). This means that the 
quark phase is favored rather than the hadron phase in 
the mixed phase for a given density at finite tempera- 
ture. Hence, the mixed phase contains more quarks at 
finite temperature. The free energy is decreased by the 
inclusion of more quarks and this change exceeds the in- 
crease of the kinetic energy As a result, EOS of the 
mixed phase becomes softer at finite temperature. We 
can see this behavior in Fig. 2J It shows the density 
profiles within 3D cell (quark droplet) for ns = 2no 
at zero temperature (left panel) and finite temperature, 
T = 30 McV (right panel). Clearly, the volume fraction 
of quark matter at finite temperature is larger than that 
at zero temperature. 

Another point is that the onset density of the mixed 
phase is shifted to a lower value (compare the left panels 
of Fig. [5] with the right panels). This behavior also comes 
from the fact that the quark phase is favored rather than 
the hadron phase for a given density at finite tempera- 
ture, as discussed above. 



mixed phase is much different from those of pure hadron 
matter. Actually hadron matter is positively charged 
and its density becomes lower in the mixed phase, both 
of which disfavor the mixture of hyperons. At finite tem- 
perature hyperon mixture becomes rather easy since the 
Fermi seas of hadrons are diffused by the thermal effect 
and entropy is increased [42| (see lower right panel); the 
hyperons appear at low density (~ no) at T = 30 MeV, 
while they do above 2uq at T — (the lower left panel). 
However, we can see that hyperons are also suppressed in 
the mixed phase even in such a case for the same reason 
as at zero temperature (see upper right panel). [55j 

C. Maximum-mass problem 

At the last of this section, we briefly discuss some im- 
plications of our results on the maximum-mass problem. 
It is well-known that EOS becomes too soft to support 
the canonical mass of ~ 1AM@, once hyperons are taken 
into account above the nuclear density [431 ]. The decon- 
finement transition to quark matter has been suggested 
to circumvent the difficulty [HI, [43| . We have confirmed 
this suggestion in the previous paper [l8| by treating the 
hadron-quark mixed phase in a proper way. This argu- 
ment has been done for T = case, but similar discussion 
is possible within the present framework. We show the 
mass - central density and mass - radius relations for 
isothermal hybrid stars in FIG. obtained by solving 
the Tolman-Oppenheimer-Volkoff equations. Below the 
sub-nuclear density (n < 0.1 fm~ 3 ), we use EOS by the 
Lattimer & Swesty [44| . We can see that the maximum 
mass M max at finite temperature is slightly smaller than 
that at zero temperature; M max = 1.42M at T = 30 
MeV and M max = 1.43M at zero temperature. This re- 
sult is easily understood by considering that hybrid stars 
in the maximum-mass region are predominantly governed 
by the EOS at high-density, where pure quark matter is 
realized; the ratio of temperature to chemical potential 
at several times the nuclear density is much less than 1 
for T < 30 MeV, which makes the EOS a little stiff. Note 
again that matter considered in this paper is neutrino- 
free, while neutrino-trapped and/or iscntropic case would 
be more interesting for, e.g., supernovae. However, these 
figures may give us some suggestions about the stability 
of newly born neutron stars after the mergers of neutron 
star binaries. 



B. Hyperon suppression in the mixed phase 

Particle fractions of quark and hadron species are 
shown in Fig.[5j upper panels show them during the phase 
transition, while the lower panels for hypothetical pure 
hadron matter for comparison. In the previous papers 
hyperon suppression has been demonstrated in the mixed 
phase at zero temperature [l8| (see upper and lower left 
panels), since the properties of hadron matter inside the 



IV. CONCLUSIONS AND DISCUSSIONS 

We have studied the hadron-quark phase transition 
at finite temperature. We have taken into account the 
finite-size effects imposing the Gibbs conditions on the 
phase equilibrium, and calculating the density profiles in 
a self-consistent manner. 

At finite temperature, EOS of the hadron-quark phase 
transition gets close to that given by the Maxwell con- 
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FIG. 3: (Color online) The Left panel shows that the droplet radius R dependence of the free energy per baryon for ub = 2no 
and different temperatures. The quark volume fraction (R/Rw) 3 is fixed to be the optimal value at T = MeV for each curve. 
The free energy is normalized by the value at R — > oo. Filled circles on each curve shows the energy minimum. The results are 
for B — 100 MeV/fm 3 , a — 40 MeV/fm 2 . The Right panel shows each contribution to the R dependence of the free energy, 
the Coulomb energy, the surface energy or the correlation energy per baryon (E c /A, E s /A, E COTT /A) at T = MeV(thick lines) 
and T = 50 MeV(thin lines). For the definite meaning of the correlation energy, see the text. 
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FIG. 4: (Color online) Density profiles and Coulomb potential c/> within 3D (quark droplet) for ns = 2no at zero temper- 
ature (left panel) and finite temperature, T — 30 MeV (right panel). The cell sizes are Rw = 23.7 fm (T = MeV) and 
Rw = 28.8 fm (T = 30 MeV). The droplet radii are R = 12.5 fm (T = MeV) and R = 18.7 fm (T = 30 MeV). 



struction. It is due to the mechanical instability of the ge- 
ometrical structure induced by the thermal effect. Pres- 
sure of the mixed phase at finite temperature are 10-30 
% smaller than that at zero-temperature th oug h the sim- 
ilar behavior appears without hyperons [HI, [33( ■ This be- 
havior is characteristic of the hadron-quark mixed phase 
which can be explained by the large degree of freedom 
in the quark phase. Hypcron fractions are suppressed by 
the appearance of the mixed phase, as in the case of zero 
temperature. 

Our calculations are subject to the neutrino-free (low 
Yi) case at finite temperature. Such situation will appear 
in mergers of neutron star-neutron star binaries or black 
hole-neutron star binaries [45| , and our result may change 
their dynamical aspects. Of course, we should take into 



account isentropic and Yi ^ situation in the core of 
supcrnovac. This work is now in progress. 



Finally we note again that EOS has many uncertain- 
ties, especially for quark matter. We simply adopted the 
thermodynamic bag model in this paper, and used the 
density-temperature independent bag constant and sur- 
face tension, while it would be interesting to include such 
dependence for a realistic description [Zfl [47j . Moreover, 
chiral restoration or color super conductivities may also 
change our results [H, [H, [43, [EO, HH • These are open 
questions for astrophysics and nuclear physics. 
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FIG. 5: (Color online) Particle fractions at zero temperature (left panels), and the finite temperature, T = 30 MeV (right 
panels). The upper two panels include quark phase and the mixed phase as well as hadron phase. Lower panels: Particle 
fractions in the hypothetical pure hadron matter. 
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